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1. Introduction 



Thermodynamical Bethe-anzatz (TBA) is one of the main methods of studying statisti- 
cal properties of integrable low-dimensional systems. Its formulation yields the functional 
equations 0, |J which have the form of multidimensional recurrent relations. The 
same equations (identical in form but different in the analytical properties of solutions) 
are known in the theory of nonlinear waves as the Hirota equations. Discrete Hirota 
equations with specified boundary conditions may yield all the known exactly integrated 
(1+1)D equations at the continuous limit. The problem of intcgrability of discrete equa- 
tions was discussed in recent papers Pi P^Oj PI. In Rcfs. ||, [|, 0, [|, [l3) a significant 
progress was achieved in applying Bethe-ansatz (BA) for integrable discrete (1+1)D sys- 
tems with the internal degrees of freedom (see also jig])- It was found out thereby that the 
discrete time variable and the variable corresponding to the internal degrees of freedom, 
as well as discrete rapidities enter the Hirota equation as the arguments of the searched 
function on equal footing. The limit of continuous time in this case corresponds to large 
spin values ]Tr| . 

The behavior of the considered system is universal due to several reasons. The for- 
mulation of the model on the lattice leads to the solutions 0, which are expressed via 
elliptic functions. Complexification of the parameters includes the possibility to consider 
the quantum-group situation with the deformation parameter being a root of unity |l7|, |l8) . 
Different values for the shifts of rapidities along the imaginary axis in the searched func- 
tions of the discrete equations lead to different functional equations appearing in the TBA 
approach. In the particular case, when the rapidities are much greater than these shifts, 
we deal with recurrent relations |l^, |2^, |2l| existing in the theory that takes into account 
the Haldane's generalized principle of exclusion statistics ]22| . 

The eigenvalues T^(u) of the transfer matrix, satisfying Eq. ((HJ) with certain boundary 
conditions (see below), are characterized by analytical properties of solutions specific for 
the problem under consideration. Thus, for example, the functional recurrent relations of 
a thermodynamical Bethe-anzatz [|J, being outwardly similar in form with equations for 
A\— case with zero boundary condition, are satisfied by functions with different analytical 
properties. The distinction in analytical properties consists in the fact that for recurrent 
TBA equations the shift of rapidities occurs along the imaginary axis. 

In the present paper we consider the case of real shifts along the axis of rapidities. 
The main goal of the paper is search for integrals of motion following from the analysis of 
the Hirota classical equation for zero and quasiperiodic boundary conditions in the case of 
Ak-i algebra. The integrals of motion are useful in seeking solutions of BA equations. As 
the example of the equations that require integrals of motion to be analyzed we consider 
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equations in the discrete and continuous limits for the k = 3 case. The general approach 
for studying the integrals of motion for discrete sine-Gordon model was used in Ref. |23| . 
The representation of the zero curvature p4], 25| for the classical discrete sine-Gordon 
model via the R- matrix was studied in Ref. p6|. 

The present paper is organized as follows. In the second section for consistency we 
provide the main expressions from Refs. [|l3|, H. Based on these equations we write 
out the discrete equations of motion for the non-Abelian case k = 3 under zero and 
quasiperiodic boundary conditions. Here we give their form in the continuous limit. The 
third section is devoted to the analysis of integrals of motion in the case of quasiperiodic 
and zero boundary conditions. In the discussion we analyze the form of solutions for the 
Ai algebra with the quasiperiodic boundary conditions. 



2. Hirota Equation 



The Hirota bilinear difference equation may be written down in different forms. One of the 
presentation of this equation for the function r(n, I, m) depending on the discrete variables 
n, I and m is 

arin, I + 1, m)r(n, I, m + 1) + Pt{ti, I, m)r(n, I + 1, m + 1)+ 

7r(n + 1 ,/ + l,m)r(n - l,l,m+ 1) = 0. (1) 

The arbitrary constants a, /?, 7 in this equation are restricted by a + (3 + 7 = 0. Further 
specification of Eq.(]l|) is performed with the aid of boundary conditions that set up the 
correspondence with the searched model in the continuous limit. Preserving integrability, 
various methods of the realization of the continuous limit permit to obtain the whole 
spectrum of soliton equations in the continuous case. The transform of the function 
r(n, I, m), 

T ( n 1 m ) = } a/ >> T a m ) (2) 

eliminates the arbitrary constants a,P,j in Eq. ([!]): 
T n (l + l,m)T n (l,m + 1) - r n (l, m)r„(/ + l,m+ 1) + r n +i(l + 1,to)t„„i(7, m— 1) = 0. (3) 
Substitution of the variables 

a = n , s = I + m , u = I — m — n , 
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r n (l,m) = T/; m (I - m - n) = T s a (u) (4) 

leads to the following bilinear functional relations (expressing the fusion rules) for the 
eigenvalues T"(u) of the transfer matrix in integrable lattice models for the A^-i algebra: 

T s > + l)2?(u - 1) - T^ +1 (u)T^_ 1 (u) = T? +1 (u)Tr\u) . (5) 

The discrete index a = 0, 1, . . .k for the A^-i algebra and the discrete time s in this 
equation are, respectively, the length and the height of the Young rectangular diagram, to 
which the eigenvalue T"(u) of the transfer matrix corresponds; u in Eq.(|^) is the discrete 
spectral parameter. Equation (g) preserves its form if the function T"(u) is multiplied by 
the product of four arbitrary functions Xi depending on one variable with the following 
arguments: 

T s( u ) Xi(« + u + s )X2(a -u+ s)x3(a +u- s)x4(a -u- s)T s Q (u) . (6) 
For the function 



Tg +1 (tt)Tg° -1 (u) 

which is invariant with respect to the transform M , Hirota equation (0) has the form 



y(« + i)y (« - 1) = ( 1 + ^ (tl) ) + (8) 

Y A U+l)Y A U L) {l + {Y a + l {u)yl){l + {Y a-l {u)yl y W 

The zero boundary conditions for the function T"(u) are as follows |l3| 

T s °(u) = for a < and a > k , (9) 
T*(u) = 4>{u -s-k), T°(u) = ct>(u + s) , 

N 

where <fr(u) — Y[ a ( r l( u ~yi))i a is the Weierstrass function, being its roots; N is the size 

i=l 

of the system and rj is a parameter. In the case of A±— algebra, these boundary conditions 
lead to the discrete Liouville equation |p7| . In terms of the gauge- invariant variables it 
has the form 

Y s a (u - l)Y/(u + 1) = (1 + Y« +1 (u)) (1 + J7_ x (u)) (10) 

with the boundary condition Ko(w) = 0. In the continuous limit after replacing Y^(u) = 
5~ 2 exp(-<f>(x, t)), n = x/S, s = t/5 with 5 -► Eq.© transforms to Liouville equation 
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The quasiperiodic boundary conditions for the function 

Tg(u) = e a \- 2a T? +2 (u + 2) (12) 

in the Ai-algebra case lead to a discrete sine-Gordon equation|l3| for the gauge-invariant 
function X^(u) = -A(l + Y s a (u)) @ |§ |o|: 

X a {u)x a (u) _ (A + *>+!)) (A + *,>-!)) 

As+lWs - lW - (1 + XX-(u + 1)) (1 + \X«(u - 1)) • {U > 

Due to the boundary condition (p"2|), the functions X" +1 (tt) with different values of 
variable " a" are related as 

X a s {u)X a s +1 {u + \) = l. (14) 

We consider discrete equations of motion for the A 2 algebra for the zero and quasiperi- 
odic boundary conditions, and their continuous limit. Applying the general approach of 
the paper JL3| in the case of zero boundary conditions for gauge- invariant functions (u) 
with a = 1, 2 we obtain the following system of discrete equations: 

(1 + Y 2 {u))Y}{u + l)Y s \u + 1) = (1 + Y' +1 (u))(l + YU{u))Y 2 {u) , (15) 
(1 + Y}{u))Y*{u + l)Y 2 (u + 1) = (1 + Y s 2 +1 (u))(l + Yl^uW^u) . (16) 

The Eqs.(15),(16) are the first result of this paper. 

Similarly to the consideration of the transition from Eq. (^) to Eq. (pi) we now intro- 
duce the functions <f>i(x,t) with I = 1,2 in such a way that Yg(u) = J exp(— </>;(x, i)). 
Let u — x/5, s — t/S. Then the continuous limit 5 — > yields 5/7(3) Toda equations plf : 

9201 ^Jte^-e* (17) 



<9t 2 dx 2 

'2 



<9 2 <fe <9 2 



<9i 2 <9x 2 



= 2e*»-e* 1 . (18) 



In the case of non-Abelian discrete sine-Gordon model the quasiperiodic boundary 
condition ( |l2| ) should be generalized for the Ak-i algebra in the following way: 

T s a (w) = e a \- ka T? +k (u + k). (19) 

Here the functions I°(w) satisfy the relation 

X*(u)X a s +1 (u + 1) . . . A s a+fe - 1 (ii + *-!) = (-l) fc , (20) 
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which generalizes Eq.(14). 

Consider the case k = 3. The Hirota equation for gauge-invariant functions X} (it) for 
the A2— algebra with boundary conditions ( |l9| ) can be written as 

(A + X^))*^)^.^) - (1 _ Ax i (u+1)X 2 (u + 2)) ( 21 ) 

(A + ^Wj^W^.W = ( 1 _ Ax i( M _ 2 )x|( M -l)) (22) 

This equations of motion is the second result of our paper. 

The continuous limit in these equations can be achieved as follows. Assume that A — ► 
and s — > 00 in such a way that s = t/vX for a fixed t. Similarly, assume that A — * and 
u — ► 00 with u = x/yX for a fixed x. Being parametrized with the aid of the function 
X l (x,t) = — exp(— (pi(x,t)) Eqs. ( pl| ) and ( p2|) transform in the continuous limit to the 
following equations of motion 

d 2 ^ &<k =2e *_ e *»_ e -«i+**) i (23) 



dt 2 dx 2 
d 2 4>2 d 2 4> 2 



dt 2 dx 2 



= 2e 02 - e* 1 - e - (01+02) , (24) 



which correspond to SL/(3)— case of affine Toda equation. 

Let us assume that in Eq.(21) Xj(u) — X 2 (u) = X s (u). Then we have 

X s+ i(n)X s _ x (u) _ (A + X s (u - 1))(A + X s (u + 1)) 



X s {u) (1 - XX s {u - 2)X s {u - 1))(A + X s {u)) 

or from Eq.(22) with the same assumption 

X s+1 {u)X s ^(u) = (A + X 3 (u-1))(A + X a (u+1)) 

X,(u) (1 - AX s (u + 2)X s (u + 1))(A + X s {u)) 



(25) 



(26) 



In this case we have two versions of discrete equation. These equations are symmetrical 
to each other with respect to changing the variable u — ► —u. One can symmetrize these 
equations adding together Eqs. (26), (27) 

X s (u) 



X s+ i(u)X s -i(u) 



6 



(A + X s (u))(l - (X/2)X s (u - 2)X s (u - 1) - (X/2)X s (u + 2)X s (u + l)) 



(X + X s (u-l))(X + X s (u+l)) 
The continuous limit of the Eq.(27) yields 

d 2 (f> d 2 (j> _ 
dt 2 dx 2 



e" 20 , (28) 



which is the Bullough - Dodd equation. 

Another discrete variant of the Bullough-Dodd equation was found from geometrical 
point of view in Refs. p2[ . Generalizing the viewpoint, we may conclude also that Eqs. 
( pl| ) and ( p2| ) are the discrete analog of the equations of motion in SU (3)-case of affine 
Toda model. To study solutions of the discrete equations of motion ([l5]) and ( |l6| ) and 
( |2l| ) ( ]22^ ) with zero ^ and quasiperiodic ( |l9| ) boundary conditions, respectively, one has 
to find the integrals of motion. We consider this problem in the next section. 



2. Integrals of motion 

Integrals of motion are known to play the key role in the study of types of motion. Un- 
der integrals of motion we assume such combination of the functions T"(u) which doesn't 
depent on argument s or u. In Ref. fl3]| integrals of motion for discrete dynamics were 
obtained for the A\ algebra and zero boundary conditions. A general determinant repre- 
sentation of integrals of motion was also written out for the Ak algebra with zero boundary 
conditions. In this section we shall show (i) how to get integrals of motion in the A2 case 
from the known integrals for the A\ algebra. In fact, we use this example to show how 
integrals of motion are transformed in the Backlund flow. Besides, (ii) we write down 
some integrals of motion for quasiperiodic boundary conditions. 

We shall use the equations of the linear problem, which ensue from the gauge and 
dual transformations on the lattice |Q of the values for the variables in question. These 
equations are equivalent to those resulting from the zero curvature condition p4| , |25| . 
Using the notations of Ref. [Q , the equations of the linear problem are written as follows 

g n (l, m)r n {l + 1, m) - g n (l + l,m)r n {l,m) = c _1 g n _i(J, m)r n +i(l + 1, m), (29) 
g n - 1 (l,m)T n (l,m+ 1) — g n -\{l,m+ l)r n (l,m) = cg n (l,m)T n -x(l,m+ 1). (30) 

Here c is the arbitrary constant. When comparing Eqs .(p9| ), ( |30| ) with the equations of the 
linear problem of the paper Jl3), the transformation (|2|) and the relation c + c_ = — 7/a, 



c should be taken into account in (E9j), (30). This relation is the necessary condition 
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for the function g n (l,m) to satisfy also Hirota equation. The inverse equations of the 
linear problem have the form 

g n {l,m)T n (l - l,m) - g n (l - l,m)T n (l,m) = -c~ x g n ^i(l - 1, to)t„ + i(Z, m), (31) 
g n (l,m)T n +i(l,m- 1) - g n {l,m- l)r n+ \{l,m) = -cg n+1 (l,m- 1)t„(Z,to). (32) 

In the terms of the functions T®(u) and G"(m), using the notations from (Q), when 
r„(Z, m) = Tu_ (I — m — n) = T s °(u) and m) = G™ +TO (' — m — n) = G"(it), one can 
write Eqs (|2^ - f52|) in the following way. For the direct equations of the linear problem we 
have 

G a s {u)T« +l {u + 1) - G a s+1 (u + 1)T S » = c-^-^u + l)T s a +l(u) , (33) 

G a s -\u + l)T? +1 (u - 1) - G a s -{{u)T:{u) = cG a s {u)T^(u) . (34) 

Similarly, the inverse equations are written as 

G a s {u)T«_ x {u - 1) - GU(u - 1)T S » = -c-'GtzHu^+Hu - 1) , (35) 

G»T;+V) - G a s _ l {u + l)T«+\u 1) = -cG;+ 1 («)2?(«) . (36) 

To make the comparison of the approaches and notations of various papers Jl3[ easier, 
we write out the equations of the linear problems having expressed the functions included 
in them via the eigenvalues of the transfer matrix, used in Ref. jl3). We note that the 
complete equivalence is achieved at c = — 1 and the substitution s — > — s. 



Let us focus our attention on the symmetry of pairs of equations (33), (34) and (|35|), 
©. It is seen that after the substitution T s a (u) -> GZ a s (-u) and G a s {u) -> Tl?(-u) Eqs. 
@, (H) transform to ©, @. 

In the case of zero boundary conditions @ for the function T s a (u), the function G°(m) 
satisfies also the same boundary conditions, but with the maximum values a less by a 
unity, i.e., k should be replaced by (k — 1). This feature expresses Backlund transform 
which realizes the transform — > A^-i- It was used in the paper fis]] to get nested 
Bethe ansatz equation. This transform for k = 3 case can be schematically presented in 
the following form: 

1 

Qi(u + s) Qi(u-s) 

(37) 

Q 2 {u + s) Gl(u) Q 2 (u-s) 

(j)(u + s) T}(u) T s 2 (u) 0(it-s) 
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Here Q a (u) are some boundary functions expressed via Weierstrass functions a(i]u)), 
Qa(u) — Q a (u — a), (f>(u) — <p(u — k). It follows also from the boundary condition 
(|) that Qk(u) = <p(u) (for details see Ref.[13]). 

The functions at each level satisfy Hirota equation. There is a relation between two 
adjacent levels, defined by the equations of the linear problem. The function Gl(u) corre- 
sponds to the case of Ai algebra. The equations of the linear problem written for it allow 
to obtain integrals of motion, which for it-dynamics |13|] are: 



h(s) 



his) 



Q 2 js - 2)G\_ u+1 ju) + Q 2 (a)GU-i(" ~ 2 ) 

Q 2 js + 2)Gj_ s+1 (u) + QMGj-s-^u + 2) 
Gl_ s (u + 1) 



(38) 



(39) 



For Ai case Eq. (p3|) of the linear problem for a = 1 and a = 2 may be written in the 
following way 



Q 2 (u + s)Tj +1 (u) - Q 2 (u + s + 2)Tl{u - 1) = c -1 Gi +1 (u)^(u + s) 



(40) 



G l s {u)T 2 s+l {u) - Gl +1 (u + 2)T^u - 1) = c" 1 ^ - a - ljT^u) . 
Similarly Eq.(|34|) is written down 

GKujT^^u) - Gj+iCu - 1)T» = cQ 2 ( M + s + l)T^u - 1) , 



(41) 



(42) 



Q 2 iu - s)T 2 s+1 iu -1)-Q 2 (u- S - 2)T s 2 (,i) = cG\ +1 iu)^u + s) . (43) 

Having expressed from ( [40| ) the necessary values G\ (u) and substituting them in (^) we 
get the first integral of motion 



7i(a) = ^(a-2) 



Q 2 (s - 2) [Tj_ M+1 (u)Q 2 ( S ) - - l)Q 2 ( a + 2)] 

<Ks) [T/_ u (u - 1)Q 2 ( S - 2) - T^iu - 2)Q 2 (s)] 



gg(f) [Tl-u-iiu ~ 2)Q 2 (5 - 4) - T}_ u _ 2 ju - 3)Q 2 (s - 2)] 
0(s - 4) [Tl_ u {u - l)Q 2 (a - 2) - (« - 2)Q 2 (s)] 



(44) 
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Similarly having expressed from ( ]4l| ) the corresponding values Gl(u) and substituting 
them in we get the second integral of motion for u— dynamics 



I 2 {s) = 4>(-s-2) 



Q 2 (2-s) \T*(u-l)Q 2 (-s) 



2)] 



\T 2 +u {u)Q 2 {2 - s) 



T 2 

± s+u- 



-i(«+l)Q 2 (-a)] 



<5a(- 



fc 2 +„-i(^ + 1)02(4 - a) - Tl +U _ 2 {u + 2)Q 2 {s - 2)] 



<K*- 4) [T S 2 + »Q 2 (2- S ) 



l)<3a(-*)] 



(45) 



In order to obtain integrals of motion for s— dynamics one should use the other pair of 
equations - Eqs.([4l|), (jf3|). These integrals of motion are related with the form (||) of 
boundary conditions. Here one of them is associated with the boundary condition at the 
left end of the segment [0, k] and the other - at the right end. The chosen example k = 3 
illustrates the method of obtaining integrals of motion from the integrals of motion for k 
less by a unity. 

Now we consider quasiperiodic boundary conditions. In order Eqs.([33|), (|34|) and j35|), 
(36), are fulfilled, the function G"(u) in the case of the quasiperiodic boudary conditions 
(19) should also satisfy the same boundary conditions. Therefore there are no Backlund 
flows in this case and the transition from algebra to A^-i algebra does not take place. 
It follows from the condition, that two functions Tfiu) and G"(u) satisfy one and the 
same Hirota equation with equal boundary conditions, that they are proportional to each 
other. The relations between the indices in these functions are defined by the law of the 
transformation of direct equations of the linear problem ( |33| ) , ( [34] ) into the inverse ones 
(p5|), (|3"6j). The coefficient of proportionality should retain Hirota equation invariant, i.e. 
it should coincide with the gauge transformation (Q). Thus we may write down 

G s( u ) = Xi(a + u + s)x 2 {a -u + s)x3(a + u- s)Xi{a -u- s)TI°(-u) . (46) 

We use Eq.pS) in (pi), (pi)- As the result we obtain 



(47) 



X4(fl U ;Tz!{-u)T» x (u + 1) - Xlia ; U + S ,\ 2) T_-:_ 1 (- U - WOO 



Xi(a-u-s-2) sy s+iv xi(a + u + s) 

iX 2 (a-u + s-2) a+ i 



Xi{a - u + s) 



-Tir\-u-l)T^(u), 



X2(a- 



Xi{a - u + s) 



-T~ 



\-u-l)T: +1 (u-l) 



X3(a + u 



X3(a + u- s) 



i] Tz:+i{-u)T:{u) 
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c X4(a _ u _ s _ 2 ) T - HO T H-i W • ( 48 ) 

We use this form of equations to get integrals of motion. Here we take into account 
that quasiperiodic boundary conditions limit the type of the functions xi an d X3- They 
should be periodic with the period 2k, while the functions X i an d Xi may be arbitrary. We 
single out the functions xi an d X3 from Eqs. (Et]), (Esh , having constructed the relations 

Xi (a + u + s + 2) X4 (a - « - s)Tr s °(-ix)T s a +1 (u + 1) 



Xi (a + u + s) X 4(a - u - s - 2)T^_ 1 (-u - l)T«(u) 
_ lX2 (a-u + s- 2)Ti: +1 {-u - l)T£l(u) 



c 



X2 {a-u + s)TI«_ x {- U - l)T«(u) 
X3 (a + u- S -2) _ X2(a-u + s- 2)T^ +1 (-u - 1)T° +1 ( M - 1) 

X3 (a + M-s) x2(o-t» + *)rr.°ji 1 (-«)2?(«) 



(49) 



x 4 (a-^-«)Tr°(-«)r°+ 1 1 ( M ) 

C X 4 (a - u - « - 2)T^+l(-u)T-(u) ' 



It follows from the periodicity of the functions Xl an d X3 that the right parts of 
Eqs. (|49|) , ( pQ ) are periodic with the period 2k with respect to three indices. After the 
replacement u — > u — s in Eq. ( f49| ) and the replacement u — > u + s in ( |50| ) the right part 
of these equations does not depend on the time s. Consequently, after this replacement of 
variables the right hand size of the equations is the integral of motion for s— dynamics. The 
presence of the arbitrary functions Xi an d the constant c in the formulas for the integrals of 
motion shows the degrees of freedom up to which we can consider the restriction imposed 
by them. 

Similarly we may obtain the preserving values of u— dynamics, substituting s — > s — u 



in Eq.(49) and s — * s + u in Eq.(50). For simplicity we assume that X2 = X4 = 1 a nd 
we consider the case of A%— algebra. For a — and a = 1 from Eq.([49|) we may get two 
invariant combinations of functions: 

Xi(u + s + 2) T° s (-u)T° +1 (u + l) ^ TU-u-VT^ju) 
Xi(u + s) IV^-u-ITOu) ° TVi(-«-l)T s °W' 1 ' 
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Xi(u + s + 3) _ Tj s (-u + 2)T} +1 (u + 1) _ 1 T° s (-u-l)T^ +1 (u-2)(Xe a ) 2 
Xl (u + s + l) T^i-u + l^iu) ° T^-u + ^T^u) 

The combinations made of the functions Xii which follow from (|50|), have the form 
X3 (u -a-2) _ T^j-u - 1)T° +1 ( M - 1) T° a (- M )r fl V 2 ( M + 2) 



Xafo -£-1) T° s (-u - l)T} +1 (u - 1) Tl a (-u + 2)T° +1 (u) 

xs(u-s + i) TViMW c rVi(-^°(u) • 1 > 

The Eqs.(49)-(54) are the main result of the paper. 

It is seen directly from these equations that after the substitution u — > u — s (u — > u+s) 
or s — > s — u (s — > s + it) in Eqs.(^l|),(|5^) ( (|53[) ,[54|) ) we obtain in the right hand sizes of 
these equations the magnitudes which do not depend on s (on u) after the change of these 
discrete variables. 



4. Discussion 

The preserving combinations (|5|) - ( |54| ) of the functions T£{u) may be used for studying 
discrete dynamics in sine-Gordon model. To find the solutions of the equation of motion 
in the discrete sine-Gordon model, Bethe-Ansatz method was used in paper g. In the 
case of Ai-algebra for the function 

N 

T»>\u) = A°>V«« J] ° - *?'H»))) > (55) 

where the index a takes the values 0,1, the function /i(s) = /io« + /Xi (/Uo>Mi are t ne 
arbitrary constants), two functions 2^-' (s) are the zero functions of Weierstrass function 
and are defined by the solutions of BA equations: 




cr(r,(zH a + l)-z^s + l)-l)) \ ( <t(t,(zI(s)-z$(s) + 1)) \ 
a( V (zj(s + 1) - z°(* + 2) - 3)) ) \o( v {z}(s + 1) - zj(s + 1) + 1)) J 
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<j( V (z}( S + 1)-z°( S + 2)-2)) 
a( v (z}( S + l)-z°( S ))) 



A A 1 



(56) 



' <j( V (zf( S + l)-z}(s + l) + 2)) \ ( a( V (z?(s) 
+ 1)-* J 1 (* + 1) + 1)) / 



1)) 



A A 1 



(57) 



In the rational limit, when a(r)u 
equations || 

(A s ) — A S+1 A S _ 1 



a{ V (z°{ S + 2)-z)( S + 2))) A° +1 Al +1 

u, the functions A s may be found from the system of 



0n2 



e<*X\Ay = a{A») 



{A° s f-A° s+1 AU=e- a \ 2 {A\Y^b{AlY 



(58) 
(59) 



The equations (56) - (|5^) describe the discrete dynamics of the zeroes Zj' (s) as the 
function of the time s. 

In the presence of the restrictions following from (p3l|)-(|5l|) to solve Eqs.(|56|), j57|), 
one should find the solutions of Eqs.(|58|), (p9[). We pay attention to the invariance of 
Eqs.(|58|), ( |59| ) with respect to the transformation of (A®,Al) — > (A°, Al)e Ss , where 5— is 
the constant. This invariance is already displayed in Eq.(^5|) with the aid of the functions 
exp(/i(s) u) with /i = ^jlq s + /m.. 

It follows from Eqs.(|58|), ( [59"1) that the functions A®' 1 have the form 



A" 



Ae~ ls 



Al =Ae^ s2+a/2 , 



(60) 
(61) 

Here the coefficients 7 and A are connected by the relation A 2 = 1 — exp(— 27), A is the 
constant. Using the solutions obtained (|60|), ( |6l| ) in the right hand size of Eqs.(|56|), d57j ) 
we get that the right hand size of these equations equals — 1. In this case, we have to 
change certainly the functions in the left hand size of Eqs.(56),(57) by a{rju) — ► u. 

Substituting the solutions r s 0,1 (u) of (55) to the gauge-invariant functions Y®(u) we 
see that taking into account ( |60| ) the dependence of these functions on s is defined only 
by the dependence of the roots z®' 1 on the discrete time s and does not depend on Gauss 
multipliers (56), (57). If we substitute the functions T s 0,1 (it) from Eq.(|55"|) to the limitations 
(|5l|)-(|5^) expressing the conservation laws in the considered system, this will allow to find 
the limitations on the dependencies of the roots z®' 1 on the discrete time s. In the case 
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N = 2 this problem may be presented in the form of the finite number of mappings. The 
thermodynamic limit N — > oo requires numerical calculations. 

In conclusion, we found the integrals of motion for s— and u— dynamics for Ak— algebra 
with quasiperiodic boundary conditions. The case k = 2 corresponds to the discrete sine- 
Gordon model. It is shown that for A2 — algebra the discrete equations of motion at zero 
boundary conditions convert in the continuous limit to the equations of motion of non- 
Abelian generalization of Liouville equation and at the quasiperiodic boundary conditions 
- to the equation of motion of the SU (3) affine Toda equation and to the equation of 
motion of the Bullough-Dodd model. 
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